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We propose a set of thermodynamic guidelines aimed at facilitating more robust screening of hydrogen stor- 
age reactions. The utility of the guidelines is illustrated by reassessing the validity of reactions recently proposed 
in the literature, and through vetting a list of more than 20 candidate reactions based on destabilized LiBH4 and 
Ca(BH4)2 borohydrides. Our analysis reveals several new reactions having both favorable thermodynamics and 
relatively high hydrogen densities (ranging from 5-9 wt.% H2 & 85-100 g H2/L), and demonstrates that chemical 
intuition alone is not sufficient to identify valid reaction pathways. 



I. INTRODUCTION 

The potential of emerging technologies such as fuel cells 
(FCs) and photovoltaics for environmentally-benign power 
generation has sparked renewed interest in the development 
of novel materials for high-density energy storage. For mobile 
applications such as in the transportation sector, the demands 
placed upon energy storage media are especially stringentfi 
as the leading candidates to replace fossil-fuel-powered inter- 
nal combustion engines (ICEs) — proton exchange membrane 
FCs and hydrogen-powered ICEs (Ha-ICEs) — rely on H2 as 
a fuel. Although H2 has about three times the energy den- 
sity of gasoline by weight, its volumetric density, even when 
pressurized to 10,000 psi, is roughly six times less than that 
of gasoline. Consequently, safe and efficient storage of H2 
has been identified^ as one of the key scientific obstacles to 
realizing a transition to Ha-powered vehicles. 

Perhaps the most promising approach to achieving the high 
H2 densities needed for mobile applications is via absorption 
in solids.^ Metal hydrides such as LaNisHg have long been 
known to reversibly store hydrogen at volumetric densities 
surpassing that of liquid H2, but their considerable weight re- 
sults in gravimetric densities that are too low for lightweight 
applications.** Accordingly, recent efforts^''''^'^'' have increas- 
ingly focused on low-Z complex hydrides, such as metal boro- 
hydrides, A/(BH4)„, where M represents a metallic cation, 
as borohydrides have the potential to store large quantities 
of hydrogen (up to 18.5 wt.% in LiBH4). Nevertheless, the 
thermodynamics of H2-desorption from known borohydrides 
are generally not compatible with the temperature-pressure 
conditions of EC operation: for example, in LiBH4 strong 
hydrogen-host bonds result in desorption temperatures in ex- 
cess of 300°Cj^ Thus the suitability of LiBH4 and other stable 
hydrides as practical H2-storage media will depend upon the 
development of effective destabilization schemes. 

Building on earlier work by Reilly and Wiswall,— Vajo et 
al}'^ recently demonstrated that LiBH4 can be destabilized by 
mixing with MgH2. In isolation, the decomposition of these 
compounds proceeds according to: 

LiBH4 ^ LiH + B + |H2, (la) 
MgH^ ^ Mg + H2, (lb) 



yielding 13.6 and 7.6 wt.% H2, respectively, at temperatures 
above 300°C. The high desorption temperatures are consistent 
with the the relatively high enthalpies of desorption: /S.H ^ 
67 (LiBH4) and -70 (MgH2) kJ/(mol H2)."''^ By mixing 
LiBH4 with MgH2, Aff for the combined reaction can be 
decreased below those of the isolated compounds due to the 
exothermic formation enthalpy of MgB2: 

LiBH4 + ^MgH2 LiH + ^MgBa + 2H2 . (2) 

That is, formation of the MgB2 product stabilizes the dehydro- 
genated state in Eq.|2] relative to that of Eq.[Tl thereby desta- 
bilizing both LiBH4 and MgH2. By adopting this strategy, 
measured isotherms for the LiBH4 + iMgH2 mixture over 
315-400°C exhibited a 25 kJ/mol H2 decrease in Ai? rela- 
tive LiBH4 alone, with an approximately tenfold increase in 
equilibrium H2 pressure In addition, the hydride mixture 
was shown to be reversible with a density of 8-10 wt.% H2. - 
Nevertheless, the extrapolated temperature of T = 225°C at 
which = 1 bar is still too high for mobile applications, 
and suggests that additional destabilization is necessary. 

The concept of thermodynamic destabilization appears to 
offer new opportunities for accessing the high H2 content 
of strongly-bound hydrides. However, the large number of 
known hydrides suggests that experimentally testing all pos- 
sible combinations of known compounds would be impracti- 
cal; thus a means for rapidly screening for high-density H2- 
storage reactions with appropriate thermodynamic si^ would 
be of great value. Towards these ends, here we employ first- 
principles calculations to identify new H2-storage reactions 
with favorable temperature-pressure characteristics based on 
destabilizing LiBH4 and Ca(BH4)2^ by mixing with selected 
metal hydrides. Our goal is to determine whether addi- 
tional destabilization of LiBH4 and Ca(BH4)2 — ^beyond that 
demonstrated" with LiBH4/MgH2 — is possible by exploiting 
the exothermic formation enthalpies of the metal borides. We 
focus specifically on thermodynamic issues since appropri- 
ate thermodynamics is a necessary condition for any viable 
storage material, and thermodynamic properties are not eas- 
ily altered. While kinetics must also be considered, catalysts 
and novel synthesis routes have been shown to be effective at 
improving reversibility and the rates of H2 uptake/release. — 
By screening through —20 distinct reactions, we identify four 
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new destabilized mixtures having favorable Gibbs free ener- 
gies of desorption in conjunction with high gravimetric (5- 
9 wt.%) and volumetric (85-100 g H2/L) storage densities. 
The predicted reactions present new avenues for experimental 
investigation, and illustrate that compounds with low gravi- 
metric densities (i.e., transition metal hydrides) may yield vi- 
able H2-storage solutions when mixed with lightweight boro- 
hydrides. An advantage of the present approach is that it relies 
only on known compounds with established synthesis routes, 
in contrast to other recent studies which have proposed H2- 
storage reactions based on materials which have yet to be 
synthesized.'^i^^-'^i^^i'^ 

An additional distinguishing feature of this study is the de- 
velopment of a set of thermodynamic guidelines aimed at fa- 
cilitating more robust predictions of hydrogen storage reac- 
tions. The guidelines are used to vet the present set of can- 
diate reactions, and to illustrate how other reactions recently 
reported in the literature'^ are thermodynamically unrealis- 
tic. In total, this exercise reveals some of the common pitfalls 
that may arise when attempting to simply "guess" at reaction 
mechanisms. 



n. METHODOLOGY 

Our first-principles calculations were performed using 
a planewave-projector augmented wave method (VASP)^"'^' 
based on the generalized gradient approximation^^ to density 
functional theory. All calculations employed a planewave cut- 
off energy of 400 eV, and k-point sampling was performed 
on a dense grid with an energy convergence of better than 
1 meV per supercell. Internal atomic positions and external 
cell shape/volume were optimized to a tolerance of better than 
0.01 eV/A. Thermodynamic functions were evaluated within 
the harmonic approximation,^^ and normal-mode vibrational 
frequencies were evaluated using the so-called direct method 
on expanded supercells.^'*-^^-^^'^^ Further information regard- 
ing the details and experimental validation of our calculations 
can be found elsewhere. ^^'^^'^^ 

Our search for high-density H2-storage reactions is based 
on a series of candidate reactions that are analogous to Eq.|2l 



y A(BH4)„ + MH^ ^ y AU„ + MBy 



3yn + X 



H2, (3) 



where v4 = Li or Ca [n = 1 (2) for Li (Ca)], M represents 
a metallic element, and the coefficients x and y are selected 
based on the stoichiometrics of known hydrides MH^ and 
borides MBj,„. To maximize gravimetric density we limit M 
to relatively light-weight elements near the top of the periodic 
table. In the case of A = Li, the enthalpy of Eq. [3]per mol H2 
can be expressed as: 



AH 



3y + X 



-AH 



(4) 

where A_ff* are the desorption (formation) enthalpies of the 
respective hydrides (borides) per mol H2 (A/). Thus AH for 
the destabilized LiBH4 reaction is simply an average of the 



hydride desorption enthalpies, less the enthalpy of boride for- 
mation. 



III. RESULTS 

Table U lists theoretical H2 densities, and calculated dehy- 
drogenation enthalpies and entropies for several potential H2- 
storage reactions. Reactions 1-22 enumerate the candidate 
new reactions, while reactions 23-27 are included in order to 
validate the accuracy of our predictions by comparing with 
experimentally-measured enthalpies''''^-^'' and previous first- 
principles results'^ (shown in parentheses). Turning first to 
the reactions from experiment (24-27), it is clear that the cal- 
culated T = 300 K enthalpies are generally in good agree- 
ment with the measured data. As mentioned above, reaction 
24 was studied by Vajo and co-workersi' (see Eq.|2]). Our cal- 
culated enthalpy of 50.4 kJ/mol H2 overestimates the experi- 
mental value by ~10 kJ/mol. However, since the experimental 
measurements were made at temperatures (T — 315^00°C) 
above the LiBH4 melting point (T„i — 268°C)<^) and our cal- 
culations are with respect to the ground state Pnma crystal 
structure,'' we expect AH""^" (Pnma) > AH'^'^'^^ (liquid) due to 
the higher enthalpy of the liquid state. 

We begin our discussion of the candidate reactions by com- 
menting on the vibrational contributions (AS'vib) of the solid 
state phases to the total dehydrogenation entropy, AS. Based 
on the notion that AS is largely due to the entropy of H2 



[AS* 



S"' 



130 J/(mol K) at 300 K], a dehydrogena- 



tion enthalpy in the approximate range of 20-50 kJ/mol H2 
would yield desorption pressures/temperatures that are con- 
sistent with the operating conditions of a FCj^ However, as 
shown in the last column of Table |I] the calculated AS'vib are 
not negligible (up to 2 1 %) in comparison to S^"^ , calling into 
question the assumption AS" ~ S^'^ and the guideline AH 
= 20-50 kJ/mol H2. This suggests that a precise determi- 
nation of the pressure-temperature characteristics of a given 
desorption reaction requires evaluating the change in Gibbs 
free energy [AG(r)], accounting explicitly for the effects of 
temperature and AS'vib, as done below. 



A. Thermodynamic Guidelines 

A key concern when attempting to predict favorable hydro- 
gen storage reactions is to ensure that the thermodynamically 
preferred reaction pathway has been identified. This is a non- 
trivial task, and our experience has shown that intuition alone 
is not sufficient to correctly identify realistic reactions involv- 
ing multicomponent systemsi^^ In this regard, several of the 
reactions in Table U (denoted by *) are noteworthy as they il- 
lustrate the difficulties that may arise when "guessing" at re- 
actions. For example, all of the candidate reactions are written 
as simple, single-step reactions. While this may seem reason- 
able given the mechanism proposed in Ref. 1 1 (Eq. |2]) and its 
generalization in Eq.[3j as we discuss below, some of these re- 
actions should proceed via multiple step pathways, with each 



3 



TABLE I: H2 densities and calculated thermodynamic quantities for candidate H2 storage reactions. Units are J/K/mol H2 for AS'vib and 
kj/mol H2 for AE and AH\ column 7 refers to the temperature at which Ph^ — 1 bar. Reactions denoted with a * will not proceed as written 
(see text). The enthalpies of reactions 24—27 have been measured in prior experiments, and are included here (in parentheses) to validate the 
accuracy of our calculations. For comparison, system-level targets for gravimetric and volumetric density are cited in the bottom row.^** 
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Reaction 


Wt.% 


Vol. density 
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A c'T=300K 


No. 




(kg H2/kg) 


(g H2/L) 






(°C) 
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12.4 
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15* 


rafRHiV 4- Mp^ TaHo + MpRi 4- ^Hi 


6.4 


79 


60.6 


41.9 


111 


—22.0 


16* 


Ca(BH4)2 + Al— > CaHo + AIB2 + 3H2 


6.3 


83 


71.7 


53.4 


200 


— 19.5 


17* 
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51.2 
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39 
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18 
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29.2 


—20 


— 15.9 


19 


Ca(BH4)2 + TiH2^ CaH2 + TiB2 + 4H2 


6.7 


106 


\1A 


1.1 




-17.7 


20 


Ca(BH4)2 + VH2^ CaH2 + VB2 + 4H2 


6.6 


108 


20.8 


3.8 




-16.2 


21* 


Ca(BH4)2 + CrH2^ CaH2 + CrB2 + 4H2 


6.5 


113 


29.9 


13.1 


-180 


-13.6 


22 


Ca(BH4)2 + Cr ^ CaH2 + CrBo + 3H2 


5.0 


86 


45.6 


27.2 


-38 


-16.4 


23 


6LiBH4 + CaH2^ CaBe + 6LiH + IOH2 


11.7 


93 


61.9 (63)" 


45.4 


146 


-22.7 


24 


2LiBH4 + MgH2-^ MgB2 + 2LiH + 4H2 


11.6 


96 


65.6 


50.4 (41)* 


186 


-21.7 


25 


2LiBH4-* 2LiH + 2B + 3 H2 


13.9 


93 


81.4 


62.8 (67)* 


322 


-27.1 


26 


LiBH4^ Li + B + 2H2 


18.5 


124 


103.5 


89.7 (95)" 


485 


-15.3 


27 


MgH2^ Mg + H2 


1.1 


109 


64.5 


62.3 


195 


1.3 
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Step having thermodynamic properties that are distinct from 
the presumed single-step pathway. 

We group the examples of how chemical intuition might 
fail into three categories, and for each class, give a general 
guideline describing the thermodynamic restriction: 

(1) Reactant mixtures involving "weakly-bound" com- 
pounds: We refer here to systems where the enthalpy to de- 
compose one (or more) of the reactant phases is less than 
the enthalpy of the proposed destabilized reaction; thus, the 
weakly-bound phase(s) will decompose before (i.e., at a tem- 
perature below that which) the destabilized reaction can pro- 
ceed. Two examples of this behavior can be found in Table U 
The first case pertains to reactions 13-16, which, based on 
their larger enthalpies relative to reaction 12, would appear to 
"stabilize" Ca(BH4)2. In reality, Ca(BH4)2 will decompose 
before (with = 1 bar at T = 88°C) any of the higher tem- 
perature reactions 13-16 will occur (T > 110°C), indicating 
that it is impossible to stabilize a reaction in this manner Ad- 
ditional examples of this scenario occur in reactions 1, 8, 17, 
and 21, which involve the metastable AIH3 and CrH2 phases. 



In the case of reaction 1, AIH3 will decompose first (yielding 
Al and IH2), followed by reaction of Al with LiBH4 (reaction 
2). The consequences of this behavior are significant, since 
although the intended reaction 1 has an enthalpy (-^40 kJ/mol 
H2) in the targeted range, in reality the reaction will consist 
of two steps, the first of which has an enthalpy below the tar- 
geted range (AIH3 decomposition), while the second (reac- 
tion 2) has an enthalpy above this range. Guideline 1: The 
enthalpy of the proposed destabilized reaction must be less 
than the decomposition enthalpies of the individual reactant 
phases. 

(2) Unstable combinations of product or reactant phases: 
Reaction 4 illustrates how the seemingly straightforward pro- 
cess of identifying stable reactant and product phases can be- 
come unexpectedly complex. Here, the starting mixture of 
LiBH4 and Mg is unstable and will undergo the exothermic 
transformation: 

2LiBH4+Mg-. ^LiBH4 + ^MgH2 + iMgB2 + iLiH, (5) 

which will consume the available Mg and form MgH2, which 
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will itself react endothermically with the remaining LiBH4 
according to reaction 24. The exothermic nature of Eq. (|5]l 
can be understood by noting that the enthalpy of reaction 4 
(46.4 kJ/mol H2) is lower than the decomposition enthalpy 
of MgH2, given by reaction 27 (62.3 kJ/mol H2). Therefore, 
the total energy can be lowered by transferring hydrogen to 
the more strongly bound MgH2 compound. Guideline 2: If 
the proposed reaction involves a reactant that can absorb hy- 
drogen ( such as an elemental metal), the formation enthalpy 
of the corresponding hydride cannot be greater in magnitude 
than the enthalpy of the destabilized reaction. 

(3) Lower-energy reaction pathways: Reaction 3, involv- 
ing a 4:1 mixture of LiBH4:MgH2, as well as the related 
reaction involving a 7:1 stoichiometry, 7LiBH4 + MgH2 — > 
MgBj + 7LiH + II.5H2, were recently suggested in Ref. TH 
which considered only a single-step mechanism resulting in 
the formation of MgB4 and MgBy, respectively. Here we 
demonstrate that these reactions will not proceed as suggested 
there due to the presence of intermediate stages with lower 
energies. In fact, both hypothetical reactions have larger en- 
thalpies [AE = 69 (4:1) and 74 (7:1) kJ/mol Hi'"*] than the 
2:1 mixture (reaction 24), suggesting that, upon increasing 
temperature, the 4:1 and 7:1 mixtures will follow a pathway 
whose initial reaction step is the 2:1 reaction (reaction 24), 
which will consume all available MgH2. Subsequent reactions 
between unreacted LiBH4 and newly-formed MgB2 will be- 
come thermodynamic ally feasible at temperatures above that 
of reaction 24, since their enthalpies exceed 50 kJ/mol H2. 
[Similar behavior is expected for reactions 9 & 10, as the 1:1 
mixture of LiBH4:Fe (reaction 9) will initially react in a 1:2 
ratio (reaction 10), which has a lower enthalpy.] Guideline 
3: In general, it is not possible to tune the thermodynamics 
of destabilized reactions by adjusting the molar fractions of 
the reactants. There is only one stoichiometry corresponding 
to a single-step reaction with the lowest possible enthalpy; 
all other stoichiometrics will release H2 in multi-step reac- 
tions, where the initial reaction is given by the lowest-enthalpy 
reaction^ 



Temperature (°C) 
300 200 100 50 -40 




1000/T(K"') 



FIG. 1: (Color online) Calculated van't Hoff plot for reactions listed 
in TableUl The region within the dashed box corresponds to desirable 
temperatures and pressures for on-board hydrogen storage: Ph^ = 1- 
700 bar, r = -40-100° C. 



therefore too weakly bound for practical, reversible on-board 
storage. However, the candidate reactions involving mixtures 
with ScH2"'^ (reactions 7^'' and 18) and Cr (reactions 11 and 
22) desorb H2 in P-T regimes that strongly intersect the win- 
dow of desirable operating conditions. These reactions have 
room- temperature enthalpies in the range of 27-33 kJ/mol 
H2, relatively high H2 densities (5-8.9 wt.% H2 and 85-100 
g H2/L), and achieve — 1 bar at moderate temperatures 
ranging from 26 and — 38°C. Thus, via a first-principles ap- 
proach of rapid screening through a large number of candidate 
reactions, and the careful use of thermodynamic considera- 
tions to eliminate unstable or multi-step reactions, we predict 
here several reactions with attributes that surpass the state-of- 
the-art for reversible, low-temperature storage materials. 



B. Destabilized Reactions 

In total, the preceding examples reveal that great care must 
be taken in predicting hydrogen storage reactions. Having 
ruled out the specious reactions, we now discuss the thermo- 
dynamics of the remaining reactions. Using the calculated 
thermodynamic data (Table|I]i as input to the van't Hoff equa- 
tion, = -foexp(— ■^), where Pq = 1 bar, Fig.[T]plots 
the equilibrium H2 desorption pressures of these reactions as 
a function of temperature.^^ Included in the plot is a rectangle 
delineating desirable temperature and pressure ranges for H2 
storage: -40-100°C, and 1-700 bar 

As expected, our van't Hoff plot confirms that the experi- 
mental reactions having large dehydrogenation enthalpies (re- 
actions 24-27) yield pressures P <C 1 bar, even at elevated 
temperatures. On the other hand, some of the candidate reac- 
tions, for example 5 and 19, readily evolve H2 at very low 
temperatures (consistent with their low enthalpies) and are 



IV. CONCLUSION 

In conclusion, using first-principles free energy calcula- 
tions we have demonstrated that further significant destabi- 
lization of the strongly-bound LiBH4 and Ca(BH4)2 borohy- 
drides is possible, and we identify several high H2-density re- 
actions having thermodynamics compatible with the operating 
conditions of mobile H2-storage applications. Unlike other 
recent predictions, the proposed reactions utilize only known 
compounds with established synthesis routes, and can there- 
fore be subjected to immediate experimental testing. In addi- 
tion, we provide guidance to subsequent efforts aimed at pre- 
dicting new H2 storage materials by illustrating common pit- 
falls that arise when attempting to "guess" at reaction mech- 
anisms, and by suggesting a set of thermodynamic guidelines 
to facilitate more robust predictions. 
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Experimental testing of hydrogen storage materials — many of 
which are air-sensitive — can be a slow, painstaking process. For 
example, an equilibrium measurement of the extent of hydro- 
gen desorption/uptake at a single temperature may require several 
months to complete in kinetically-hindered materials.'"' In con- 
trast, the first-principles thermodynamic calculations presented 
here (encompassing more than 20 unique reactions) were com- 
pleted in 2-3 weeks. 

This discussion assumes that the entropies of all competing re- 
action pathways are similar. Our results in Table I show that 
this is generally not the case; generalization of the above guide- 
lines to the free energies is straightforward and will be presented 
elsewhere^' . 

We neglect the LiBH4 structural transition at Ts ~ 108° G^, which 
should reduce the slope of the data in Fig.[T]for T > Ta. 
^* It should be noted that the high cost of Sc may preclude its use in 
practical applications. 
This reaction was also reported in Ref.ls^ 



